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INTRODUCTION 


Synchrotron  radiation  is  the  major  deceleration  mechanism  for 
relativistic  electrons  trapped  in  Jupiter's  magnetic  field  [ Birmingham 
et  al.  , 1974-  Coroniti,  1974],  Since  the  deceleration  time  greatly 

exceeds  the  bounce  time  for  such  electrons,  it  is  of  great  interest 
to  have  available  some  simple  (closed-form)  expressions  for  the  cor- 
responding bounce-averaged  transport  coefficients  that  appear  in  the 
Fokker- Planck  equation.  Previous  investigators  have  considered  only 
those  electrons  that  mirror  at  or  near  the  magnetic  equator  in  order 
to  simplify  the  averaging  procedure.  The  purpose  of  the  present 
work  is  to  extend  the  averages  to  arbitrary  equatorial  pitch  angles, 
i.  e.  , to  arbitrary  mirror  latitudes. 

The  present  method  of  averaging  makes  use  of  a recent  and  ex- 
cellent approximation  [ Davidson,  1976]  for  the  function 

T(y)  = ( 1 /4La)  (j)  [ 1 - y2(B/B0)]‘1/2  ds 

3 /4 

a T(0 ) - [ T(0)  - T(  1 )]y 

* 1.380173  - 0.639693y3/4  (1) 

to  which  the  particle's  full  bounce  period  2tt/D^  s (4La/pc)T(y)  is 

directly  proportional.  The  equation  of  a dipolar  field  line,  identified 

2 

by  the  dimensionless  label  L,  is  r = La  sin  0,  where  r is  the  ra- 
dial coordinate,  a^  'is  the  planetary  radius,  and  0 is  the  magnetic 
colatitude.  The  coordinate  s measures  arc  length  along  the  field 
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line.  The  magnetic  field  varies  in  magnitude  from  Bq  at  the 

equator  of  the  field  line  to  infinity  at  r = 0.  However,  a particle 

of  equatorial  pitch  angle  s sin_1y  must  mirror  at  field  intensity 

B = B^/v  and  thus  remain  trapped  between  the  latitudes  at  which 
m 0 

B = B . For  a particle  of  local  velocity  pc  (where  c is  the  speed 
m 

of  light)  the  full  bounce  period  (interval  between  consecutive  visits 

to  the  same  mirror  latitude)  is  as  given  above. 

The  rate  of  change  of  kinetic  energy  during  synchrotron  radia- 

2 

tion  can  be  written  as  Ymgc  > where  y is  the  ratio  of  relativistic 
mass  (m)  to  rest  mass  (mQ).  The  rate  of  change  of  the  sine  of 
the  equatorial  pitch  angle  is  y,  and  its  bounce  average  is  denoted 
<y>.  More  generally,  the  bounce  average  of  any  quantity  Q is 
given  by 

<Q>  = /q[1  - y2(B/B0)]'1/2  ds 

* jf  [ 1 - y2(B/B0)]"1/2  ds.  (2) 

Although  and  <y>  suffice  to  describe  the  synchrotron  "trans- 

port" of  radiation-belt  electrons,  it  also  proves  useful  (for  rea- 
sons having  to  do  with  the  canonical  Hamiltonian  formalism)  to  have 
available  expressions  for  the  bounce-averaged  time  derivatives  of 
the  first  two  adiabatic  invariants  M and  J. 

It  develops  (see  below)  that  all  four  of  the  desired  bounce 

averages  (namely  (yX  (MX  and  (j))  can  be  expressed  in 

2 3 

terms  of  the  bounce  averages  ((B/Bq)  ) and  ((B/Bq)  X Moreover, 
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the  bounce  average  of  (B/BQ)  is  analytically  related  to  T(y),  i.  e. , 
can  be  derived  from  T(y)  by  taking  integrals  or  derivatives,  if  n 
is  a positive  integer.  Thus,  the  excellent  approximation  for  T(y) 
given  by  Davidson  [1976]  can  be  manipulated  to  yield  equally  ex- 
cellent approximations  for  the  Fokker-Planck  transport  coefficients 
associated  with  synchrotron  loss  in  a dipole  field.  This  is  the  ob- 
jective that  is  achieved  in  the  present  work. 


SYNCHROTRON  LOSS 


It  is  well  known  [e.  g.  , Jackson,  1962]  that  the  rate  of  energy 
loss  by  a particle  of  charge  q is  given  by 

ym  c2  = - (2q2y6/3c)[((3)2  - (pxp)2]  (3) 

in  the  presence  of  forces  that  tend  to  change  its  velocity  pc.  The 
force  relevant  for  synchrotron  radiation  is  that  which  yields  (3  = 

IVV\ 

where  ^£2 ^ = - qB/yrngC  is  2ir  times  the  gyrofr equency. 
Thus,  it  follows  f Coroniti,  1974]  from  (3)  that 


V = (2q4  B2/3m^)(Y2  - l)sin2o 

= - (2q4B2/3m03c5)(Y2  - 1)(B/B0)3y2,  (4) 


2 2 

where  a is  the  local  pitch  angle,  i.  e.  , where  sin  a = y (B/Bq)  in 
the  course  of  adiabatic  charged-particle  motion.  It  follows  at  once 
from  (4)  that 
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<Y>  = - (2q4B02/3m03c5)(Y2  - l)<(B/B0)J>y^> 


3v  2 


where  <(B/Bq)  > is  the  bounce  average  (to  be  evaluated  below)  of 

3 

the  quantity  (B/Bq)  . 

The  change  in  energy  described  by  (4)  is  accompanied  by  a 
deflection  of  the  local  pitch  angle  a.  In  order  to  evaluate  this,  it 
proves  convenient  to  decompose  the  particle  momentum  p s ym-.pc 
into  components  p = p • B = p cos  a and  p = I p X B I = p sin  a. 

It  is  well  known  that 


1 + (p/mQc)2  = 1 + (Pn  /mQc )2  + (pi/mQc)2.  (6) 


However,  the  operative  contributions  of  p^  and  p^  to  y remain  to 
be  determined.  It  seems  evident  that  synchrotron  radiation  must 
leave  pc  cos  a unaltered,  since  a Lorentz  transformation  at  this 
velocity  along  B^  will  bring  the  observer  into  a frame  in  which  the 
electron  gyrates  about  the  field  in  a collapsing  circle  (rather  than  in 
a collapsing  helix).  The  synchrotron  radiation  in  this  frame  is  sym- 
metrical with  respect  to  the  plane  of  gyration,  and  so  there  is  no 
radiation  force  that  would  impel  the  particle  to  acquire  a velocity 
component  parallel  to  B^  in  this  frame.  It  follows  from  this  argu- 
ment that  p(  = ym^pc  cos  a in  the  original  frame  of  reference. 

The  foregoing  conclusion  about  p^  is  contrary  to  that  assumed 
by  Coroniti  [1974],  who  took  p^  = 0 without  detailed  justification. 
Thus,  the  present  decomposition  of  (6)  to  yield 


22.,.  n2  2 i 

mQ  c v y[  1 - P cos  a ] 
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r 


2 2 

differs  from  his  by  the  factor  [ 1 - p cos  a].  The  interpretation 
of  the  present  result  is  that  the  synchrotron  radiation  from  a rela- 
tivistic electron  is  somewhat  biased  toward  the  direction  of  p , 

II 

and  that  the  resulting  reaction  force  (in  the  original  frame  of  ref- 
erence) is  iust  sufficient  to  keep  constant  the  component  of  particle 
velocity  (not  momentum)  parallel  to  B. 

By  combining  (7)  with  the  foregoing  relationships  p^  = ym^pc 
cos  a and  y = ( Bq  /B ) ^ ^2  si  no  = ( Bq  /B)^2(p^/p),  one  can  derive 
the  expression 


y = (Bq/B)1  /2  | <px  /p)  - (px  /p2  ) p 1 

= (v/yy)(Y2  - 1 >" 1 [ 1 - y2(B/BQ)l  (BQ/B)  (8) 


for  the  rate  of  change  of  sin  a q.  Thus,  it  follows  from  (4)  that 


<y>  = - <2  q4y  BQ2/3  ym03c5) 

X [<(B/B0)2>  - y2<(B/B0)3>] 


in  the  course  of  synchrotron  radiation. 

2 

Since  M = p^  /Zm^B,  it  follows  from  (7)  and  (4)  that  the 
bounce-averaged  time  derivative  of  the  first  adiabatic  invariant  is 
given  directly  by 


<M>  = <(p^p^/m0B)>  = < (y  it>qC2  / y B )( 1 + y2p2sin2o)> 


- (4Mq4B2/3ym03c5)K(B/B())'i>  + (2MBQ /mQc2  )<(  B/BQ ) S>  1 , 


1 

1 


2 -2 

where  1 - (3  = y as  usual.  The  same  result  can  be  obtained 

2 2 2 

from  (5)  and  (9)  by  setting  M = (y  - l)(mQc  y /2  Bfl)  and  applying 
the  chain  rule  of  differential  calculus.  This  latter  (indirect)  method 
enables  one  to  obtain  from  the  expression  J = 2LapY(y)  for 

the  second  adiabatic  invariant  in  a dipole  field,  where  [Schulz,  1971; 
Davidson,  1976] 

Y(y)  = 2y/’1(y')"2  T(y')  dy' 

Jy 

» 2T(0)  + [6T(0)  - 8T(  1 )]  y - 8[T(0)  - T(l)]y3/4 
» 2.760346  + 2.357194  y - 5. 1 17544  y3 /4.  (11) 

It  follows  from  the  first  line  of  (11)  that  y Y'(y)  = Y(y)  - 2 T(y). 
Thus,  the  second  adiabatic  invariant  has  a bounce-averaged  time 
derivative  given  by 

<j>  = 2La{(m0c/p)<v>Y(y)  - (p/y)<y>[  Y(y)  - 2T(y)]| 

= - 4Lap(q4B02/3vm03c5){[Y(y)  - 2T(y)  ]<(B/BQ)2> 

+ y2[(\2-  l)Y(y)  + 2T(y)]«B/B0)3>|  (12) 

in  closed  form.  It  remains  to  be  shown  that  the  bounce  averages 
2 3 

<(B/Bq)  y and  <(B/Bq)  > can  be  derived  analytically  from  (1). 
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BOUNCE  AVERAGES 


It  proves  convenient  to  introduce  the  auxiliary  variable  u - y . 
Bounce  averages  of  the  required  sort  can  then  be  obtained  from  the 
realization  that 

r2(B/BQ)]1/2  ds 

= - - y2(B/BQ)]'1/2  ds 

= - i<(B/B0)>[4LaT(y)],  U3a> 

y2(  B/Bq  ) ]3  72  ds 

= |<(B/B0)2>[4LaT(y)].  (13b) 

y2(  B/Bq  ) ]5  /2  ds 

= - i|-<(B/B0)3>[4LaT(y)],  (13c) 

and  (more  generally) 

- y2(B/B0)l<2”->'/2  d, 

dun-' 

= ( - l/2)n(2n  - 1 ).'  ! <(B/BQ)n>[  4LaT(y)].  (13d) 


Thus,  the  problem  at  hand  will  be  solved  by  the  discovery  of  a 
prescription  for  evaluating  the  integral  on  the  left-hand  side  of  (13d). 


The  prescription  follows  from  the  fact  that 


£{y‘-2n/[>  - y2(B/B0!<2"-‘>'2  d.} 

= - (B/B0)]|Zn-l)/2  ds 

- - -Hsr/t1  - y2(B/B0)]<2n-3l/2  d,.  ,14) 

2y  J 

Since  the  integrals  found  on  the  left-hand  side  of  (13d)  vanish  at 
y = 1 for  positive  integers  n,  it  follows  from  (14)  that 

j>[  1 - y2(B/B0)](2n_1)/2  ds 

= y211'1/1  ^n;^2~n  (y')2(B/B0)](2n~3)/2  ds  dy' 

(15) 

for  nil.  Thus,  it  follows  from  (1)  that 

jffl  - y2(B/BQ)]1/2  ds  = y/’1(y')'2[4LaT(y')]  dy' 

= 2 La  Y(y)  (16a) 

for  n = 1,  from  (11)  that 


<f\  1 - y2(B/BQ)]3/2  ds  = 3 y2  Z*1  (y ')~4  [ 2La  Y(y ')  ] dy' 

J Jy 

« (2La/3 ) { 6T(0)  + 9[3T(0)  - 4T(l)]y 


for  n = 2,  and  from  (16b)  that 


j>  f 1 - y2(B/B0)]5/2  ds 

« 4La  | T(0)  + (15/8)  [ 3T(0)  - 4T(l)]y 

- (5/12)  [ T(0)  - 4T(l)]y3 

- (320 /5 1 ) [ T(0)  - T(  1 ) ] y3 /4 

+ (3/\3b)[3T(0)  - 20T(l)]y5}  (16c) 

for  n = 3.  It  follows  from  (13)  and  (16),  upon  insertion  of  the 

2 

auxiliary  variable  u = y , that 

<(B/BQ)>T(y)  a 3 [ T(0)  - T(  1 ) ] y"5/4  - [3T(0)  - 4T(l)]y'1, 

(17a) 

6<(B/BQ)2>T(y)  a 

10  [ T(0)  - T(l)]y_13/4  - 3 [ 3T(0 ) - 4T(l)]y’3 

+ [ 4T(  1 ) - T(0 ) ] y” 1 , (17b) 

and 

408  <(B  /BQ)3>  T(y)  a 

520  [ T(0 ) - T(l)]y”21/4  - 153[3T(0)  - 4T(l)]y"5 

+ 34  [ 4T(1 ) - T(0)]y'3  + 9[20T(1)  - 3T(0)]y-1.  (17c) 

Evaluation  of  the  relevant  transport  coefficients  is  achieved  by 
substituting  (17)  in  (5),  (9),  (10),  and  (12).  One  obtains,  for 
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example 


<Y>  a - (q4B2/6l2m3c5)(y2  - 1){  520[T(0)  - T(l)]y~13/4 

- 153  [ 3T(0)  - 4T(  1 )]  y-3  + 34[4T(1)  - T(0)]y_1 

+ 9 [ 20T(  1 ) - 3 T( 0 ) ] y | 4 T(y)  (18) 

upon  substitution  of  (17c)  in  (5).  Since  T(0)  ~ 2T(1),  all  of  the 

square-bracketed  coefficients  in  (18)  are  positive.  The  correspond- 
ing expression  for  <y>  is  given  by 

<y>  a - (q4B02/612Ym03c5)  { 160[T(0)  - T(l)]y"9/4 

- 51  [ 3T(0)  - 4T(  1 ) ] y-2  + 34[4T(1)  - T(0)] 

- 9 [ 20T(  1 ) - 3 T(0 )J  y2  | 4 T(y)  (19) 

upon  substitution  of  (17b)  and  (17c)  in  (9).  It  is  easy  to  verify  by 

direct  evaluation  of  (19)  and  (18)  that  <y^  0 and  that  <y^  ap- 

proaches the  equatorial  value  of  (4)  as  y — * 1.  Tnese  desired 

3 2 

results  are  assured  by  the  fact  that  <(B/Bq)  > — • ((B/Bq)  ^ -• 
^(B/Bq)>  — 1 when  this  limit  is  taken  in  (17).  The  bounce  average 
of  B/Bq,  as  given  by  (17a),  is  not  actually  utilized  in  the  theory 
of  synchrotron  loss  but  is  included  here  in  the  interest  of  algebraic 
completeness. 

It  is  convenient  that  the  y dependence  and  the  y dependence 
appear  in  separate  factors  in  both  (18)  and  (19).  This  factoriza- 
tion makes  the  quantities  ( - 3 m^c3  /2  q4B^)(  mQc/p)2<Y>  and 
( - 3m^ c^/2  q4Bg  )y(y)  universal  functions  of  y.  Both  are  plotted 

in  Figure  1,  the  former  as  a solid  curve  and  the  latter  as  a dashed 

2 3 

curve.  Thus,  the  solid  curve  represents  y <(B/Bq)  > and  the 
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DIMENSIONLESS  (NORMALIZED!  LOSS  RATES 
(-3rmjc5/2q4B20Ky) (-3mJ  c5/2q4  B 


EQUATORIAL  PITCH  ANGLE  aoUsirr’y),  deg 


Figure  i.  Maximal  values  ("data"  points)  and  bounce  averages  (dashed 
and  solid  curves,  respectively)  of  the  quantities  y(B/B(j)^ 

- y3(B/Bo)^  and  y2(B/Bo)^«  which  are  proportional  (respec- 
tively) to  the  loss  rates  y and  y due  to  synchrotron  radiation 
by  particles  trapped  in  a dipolar  magnetic  field 


dashed  curve  represents  y ((B/Bq)2^  * y3((B/Bq)3X  The  isolated 

2 3 

"data"  points  represent  instantaneous  values  of  y (B/BQ)  and 

O O 1 

y (B/Bq j - y (B/Bq)  , i.  e.  , correspond  to  instantaneous  values 
of  v and  y plotted  on  the  same  scale.  The  filled  circles  repre- 
sent y”4  and  correspond  to  mirror-point  values  of  y.  The  open 
circles  represent  y - y3  and  correspond  to  equatorial  values  of  y. 

The  instantaneous  value  of  y vanishes  by  virtue  of  (8)  at  the  par- 

2 3 

tide's  mirror  points.  Since  y,  being  proportional  to  y (B/BQ)  , 
is  a monotonic  function  of  B/Bq,  it  is  not  surprising  that  the 
mirror-point  values  (filled  circles)  of  |y|  exceed  the  bounce  aver- 
age (solid  curve).  However,  it  seems  that  |y|  is  a monotonic  func- 
tion of  B/Bq,  and  therefore  maximal  at  the  equator,  only  for  y2 
> 2/3.  Otherwise,  the  maximal  instantaneous  value  of  |y|  is  at- 
tained where  B = 2BA/3y2  = 2B  /3.  The  crosses  (X)  in  Figure  1 
thus  correspond  to  the  maximal  instantaneous  values  of  |y|  and 
impose  (for  y £ 2/3)  an  upper  bound  of  (4/27y  ) on  the  quantity 
y <(B/Bq)2>  - y3<(B/BQ)3>  to  which  y is  directly  proportional,  ac- 
cording to  (8)  and  (4).  Similarly,  the  open  circles  (corresponding 

-a  2 3 3 

to  y - y ) impose  an  upper  bound  on  y(( B/Bq)  y ((B/Bq)  ^ for 

y2  2 2/3. 

The  similarities  between  the  coefficients  in  (18)  and  those  in 
(19)  are  not  accidental.  There  is  a very  definite  analytical  rela- 
tionship [derivable  from  (1),  (2),  (5),  and  (9)]  between  <y>  and 
(y>.  This  relationship  follows  from  the  property  that 

( 1 /y)<y >T(y)  = - (Y2  - 1 ) (q4BQ2/12 m^La) 

X y <j)  (B/Bq)3  [ 1 - y2(B/B0)]‘1/2  ds,  (20a) 
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whereas 


( 1 /y  )<Y  >T(y)  = - (l/Y)(q4B02/I2m03c5La) 

X^(B/B0)2[1  - y2(B/B0)]1/2  ds.  (20b) 

Inspection  now  reveals  that  the  right-hand  side  of  (20a)  is  precisely 
2 

equal  to  -y(y  - 1)  times  the  derivative  (with  respect  to  y)  of  the 
right-hand  side  of  (20b).  Thus,  it  follows  from  (20)  that 


<Y>T(y) 

y ( v 2 - i)y 


d T(y) 

<y> 

dy  y 


(21a) 


and  from  (21a)  that 


<y> 


1T(y')<y(y')>  dy' 
y' (Y2  - 1)Y 


(21b) 


since  <y>  = 0 by  virtue  of  (8)  and  (9)  for  y = 1.  Inspection  of 
(18)  and  (19)  reveals  that  < y>  and  <y>,  as  calculated  via  (17), 
satisfy  precisely  the  relationship  imposed  by  (21a).  This  fact  tends 
to  confirm  the  correctness  of  the  algebraic  manipulations  carried 
out  above. 


FOKKER -PLANCK  EQUATION 

If  processes  other  than  synchrotron  loss  are  neglected  as  a 
first  approximation,  the  evolution  of  the  bounce-averaged  canonical 
phase-space  density  f is  given  f Haerendel,  1968;  Schulz  and  Lan- 


zerotti,  1974]  by 


# + wM7],  + ^[<}>7]m  = °-  1221 

Implementation  of  (22),  however,  would  require  one  to  express 
and  as  functions  of  M and  J,  whereas  the  present  work  yields 

<M>  and  <j>  mainly  as  functions  of  V and  y.  It  is  certainly 
possible  to  express  y and  y as  functions  of  M and  J [Chen  and 
Stern,  1975]  by  means  of  accurate  analytical  approximations  unre- 
lated to  (1).  However,  this  would  disturb  the  internal  consistency 
heretofore  maintained  in  the  present  work,  in  which  all  analytical 
approximations  have  been  derived  from  (1). 

The  transformation  of  (22)  to  the  variables  y and  y,  more- 
over, would  make  unnecessary  any  such  appeal  to  analytical  ap- 
proximations not  derivable  from  (1).  All  that  one  requires  in  this 

case  is  the  Jacobian  of  the  transformation  from  (M,  J)  to  (y,y). 

2 2 2 

Since  it  has  been  shown  above  that  M = (y  - lXm^c  y /2Bq)  and 
J = 2LapY(y),  where  yY'(y)  = Y(y)  - 2T(y),  it  follows  that  the 
relevant  Jacobian  is  given  by 

9(M,  J)/3(y,y)  = - 4Layp(m0c2/BQ)y  T(y).  (23) 

Thus,  it  follows  [ Haerendel,  1968;  Schulz  and  Lanzerotti,  1974] 
from  (22)  and  (23)  that 

£ + ^^[vp<7>T]y  + 97f?7^[yT,’'Ki'>T]v  * °- 
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(24) 


V 


2 1/2 

where  p = (y  - 1)  The  convenience  of  (241  is  enhanced 

by  the  fact  that  the  y dependence  and  the  y dependence  of 
and  <y>  appear  in  separate  factors  in  (18)  and  (19).  The  unit  Ja- 
cobian in  (22)  arises  from  the  fact  that  M and  J are  canonical 
action  variables  of  the  underlying  phase  space  f Haerendel,  1968; 
Schulz  and  Lanzerotti,  1974],  However,  the  forms  of  and 

O)  given  by  (10)  and  (12)  are  less  convenient  for  use  in  (22)  than 
are  the  forms  of  <y>  and  <y>  given  by  (18)  and  (19)  for  use  in  (24). 

In  the  likely  event  that  additional  dynamical  processes  (besides 

synchrotron  loss)  are  operative,  the  corresponding  Fokker-Planck 

terms  should  be  added  to  (24).  These  might  include  collisional 

terms  leading  to  energy  loss  (additional  <Y>),  pitch-angle  diffusion 

(D^y),  and  energy  diffusion  (D^ ).  They  might  also  include  wave- 

particle  interaction  terms  requiring  the  specification  not  only  of  D 

and  but  also  of  the  off-diagonal  term  D (=  D^).  Finally, 

they  might  include  distributed  sources  (S)  and  (charge-exchange) 

sinks  (-<r”^>f)  of  trapped  radiation,  as  well  as  radial  diffusion 

(D.  . ) at  constant  M and  J.  There  is  no  contradiction  in  describ- 
L L 

ing  radial  diffusion  as  occuring  at  constant  M and  J,  even  when 
the  terms  describing  other  dynamical  processes  in  the  same  Fokker- 
Planck  equation  are  being  evaluated  at  constant  y or  y.  The  same 
phase-space  density  f is  simultaneously  a function  of  (M,  J,  L;t)  and 
of  (y,y,L;  t).  The  representation  chosen  for  the  purpose  of  evalu- 
ating a given  term  in  the  Fokker-Planck  equation  should  be  deter- 
mined by  the  set  of  variables  with  respect  to  which  that  individual 
term  has  been  explicitly  specified. 
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DISCUSSION 


The  present  work  has  yielded  expressions  (in  closed  form)  for 
the  bounce-averaged  time  derivatives  of  a particle's  kinetic  energy 
and  equatorial  pitch  angle  (more  precisely,  the  sine  thereof).  The 
expressions  thus  derived  are  contingent  on  an  accurate  representa- 
tion for  T(y),  as  given  in  (1).  Experience  has  shown  [Schulz, 

1971]  that  functions  [such  as  Y(y)  in  (11)]  which  are  derived  from 
T(y)  by  integration  are  about  as  accurate  (by  percentage)  as  T(y) 
itself.  The  present  representation  [David son,  1976]  for  T(y)  is 
accurate  within  0.6 % for  0 £ y £ 1.  Since  the  bounce  average  of 
(B/BQ)n  is  found  to  involve  n integrations  over  T(y)  and  n differ- 
entiations of  the  result,  the  accumulation  of  additional  error  is  pos- 
sible  but  not  likely.  Since  the  bounce  averages  <(B/BQ)  > and 
<(B/Bq)3>  are  the  essential  ingredients  in  <Y>  and  <yX  it  is 
likely  that  <y>  and  <y>  are  thereby  represented  about  as  accurate- 
ly as  T(y)  itself,  i.  e. , within  ~ 0.6%  for  0 £ y £ 1. 

Moreover,  the  present  expressions  for  <y>  and  <y>  are  com- 
paratively simple  and  easy  to  program.  Each  is  a ratio  of  func- 
tions of  y1//4,  having  a numerator  of  four  terms  and  a denominator 
of  two  terms.  The  superficially  cumbersome  numerical  coefficients 
involve  only  T(0)  and  T(l)  in  linear  combination,  where  T(0)  « 
1.3801730  and  T(l)  « 0.7404805.  Thus,  the  expressions  for  <y> 
and  <y>  given  respectively  by  (18)  and  (19)  should  be  fully  satis- 
factory for  most  analytical  purposes. 

The  component  of  acceleration  (|3c)  associated  with  curvature 
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of  the  field  line  ( 9JB / 9s)  has  been  neglected  here,  as  in  Coroniti 
[1974],  in  its  contribution  to  y and  y.  This  approximation  should 
be  valid  for  particles  having  j (pc  /qB^ ) ( 9B/  9s  )p  | « y,  i.  e.  , for 
particles  to  which  the  laws  of  adiabatic  motion  (also  tacitly  assumed 
here)  apply. 

The  present  method  for  calculating  bounce  averages  of  the 
form  <(B/Bg)n>  should  be  useful  in  many  applications  besides  the 
theory  of  synchrotron  loss  in  a dipole  field.  The  method  lends  it- 
self readily,  for  example,  to  the  calculation  of  collisional  transport 
coefficients  in  the  presence  of  a scattering  medium  distributed  with 
a density  proportional  to  (B/Bq)11  along  a field  line.  Although  the 
method  presently  requires  n to  be  a positive  integer,  it  should  be 
easy  to  find  an  analogous  method  for  handling  negative  integers  n. 
Finally,  field  geometries  other  than  the  dipole  could  be  treated  by 
introducing  the  forms  of  T(y)  appropriate  to  the  other  field  geome- 
tries of  interest. 
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LABORATORY  OPERATIONS 


The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 

fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics:  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging:  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy:  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems, 
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